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ABSTRACT 

Low-luminosity active galactic nuclei (LLAGNs) represent the bulk of the AGN population in the 
present-day universe and they trace the low-level accreting supermassive black holes. In order to 
probe the accretion and jet physical properties in LLAGNs as a class, we model the broadband radio 
to X-rays spectral energy distributions (SEDs) of 21 LLAGNs in low-ionization nuclear emission-line 
regions (LINERs) with a coupled accretion-jet model. The accretion flow is modeled as an inner ADAF 
outside of which there is a truncated standard thin disk. We find that the radio emission is severely 
underpredicted by ADAF models and is explained by the relativistic jet. The origin of the X-ray 
radiation in most sources can be explained by three distinct scenarios: the X-rays can be dominated 
by emission from the ADAF, or the jet, or the X-rays can arise from a jet- ADAF combination in which 
both components contribute to the emission with similar importance. For 3 objects both the jet and 
ADAF fit equally well the X-ray spectrum and can be the dominant source of X-rays whereas for 11 
LLAGNs a jet-dominated model accounts better than the ADAF-dominated model for the data. The 
individual and average SED models that we computed can be useful for different studies of the nuclear 
emission of LLAGNs. From the model fits, we estimate important parameters of the central engine 
powering LLAGNs in LINERs, such as the mass accretion rate and the mass-loss rate in the jet and 
the jet power - relevant for studies of the kinetic feedback from jets. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — galaxies: nuclei 
— galaxies: jets — galaxies: Seyfert 



1. INTRODUCTION 

One of the central paradigms in extragalactic astron- 
omy is that today's galaxies host supermassive black 
holes (SMBHs) at th eir centers and have a symbiotic 



evolution wit h them (Magorrian et al. 1998 Ferrarese 



& Ford 12005). The central black holes seeded m proto 



galaxies at high redshift grow in mass during cosmic 
history through a sequence of me rgers of massive black 
holes and accretion ep isodes (e.g., [Volonteri et al 



2005 



Fanidakis et al. 2011). The massive black holes grow m 



such a way that they co-evolve with the galactic bulges 
JKormendy fe Richston e 1995; F errarese fe Merritt|2000| 
|Gebhardt et al.|200U[ |Ferrarese fe Fo rd 2005 ) but do not 
seem to d o so with the galaxy disks ( Kormen dy fe Ben- . 
der || 201 l | or dark matter haloes ( |Kormendy, Bender fe 
Cornell|[20TT|) . 

The growth of SMBHs i s dominated by the quasar 
phase ([Marconi et al.||2004[ |McLure fe Dunlop||2004| [ST 



jacki ef~al.||2007[ jMerloni fe He inz 2008) during which 
they accrete from their large gas reservoir via geometri- 
cally thin accretion disks which are radiatively effici ent 
( jShakura fe Sunyaev|[l973l |Novikov fe Thorne||1973 [) as 
suggested by different lines ot eviden ce ( Shicld sfl978| Ko- 
ratkar & Blaes 1999 Shang et al. 2005p . The quasar 
population has a strong cosmologicai evolution with its 
density decreasing by a factor of ~ 100 from z <~ 2 - 
where the density reaches its peak - to the present time 
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(e.g., |Osmer||2004| |Silverman et al. |2005[ ) . 

In the present-day universe, SMBHs are underfed com- 
pared to the ones at high z and are "sleeping". Most 
of SMBH activity at low z is dominated by the weak 
end of the AGN luminosity f unction in the form of low- 
luminosity AGNs (LLAGNs; [Ho, Filippenko fe Sargent 



1995 1997b 



Ihough the 



Nagar, Falcke fe~ Wilson||2U05| |Ho||20U. 

LAG1N phase dominates the time evoiu 



tion of SMBHs, it contributes little to their mass growth 

Hopkins fe Hernquistl | 2006l iSijacki et al. 1120071 iMerlonil 
t)l I v. ' ii, / Y„„||jvW/iI\ 



Hemz 2008; Xu fe Cao[2mTJ| ~ 
The bulk of the LLAGN population (? 



2/3; Ho 2008 



20091 reside in low-ionization nuclear emission-line re 



gent|1997al ). 

terns which 



-rr\ gions (LIN ERs; Heckman 1980 Ho, Filippenko & Sar 



LLAGNs are extremely sub-Eddington sys- 
are many orders of magnitude less lumi- 
nous than quasars, with average bolometric luminosities 
(-^boi) ~ 10 40 — 10 41 erg s - 1 and an average Eddington 
ratio of L bol /L E dd ~ 10" 5 (Ho 2009 ) where -^Edd is the 
Eddington luminosity. The observational properties of 
LLAGNs are quite different from those of more luminous 
AGNs. Regarding the broadband spectral energy distri- 
butions (SEDs), LLAGNs seem not to have the thermal 
continuum prominence in the ultraviolet (UV) - the "big 
blue bump" - which is one of the signatures of the pres- 
ence of an opt i cally thick, geometrically thin accretion 
disk jHo| |1999l |Nemmen et al.||2006| |Wu, Yuan & Cao 



2007; 



Ho pODH} |Eracleous, Hwang fe Flohic||2010[ ). Re 

garding trie emission-lines, L LAGINs typically have weak 
and narrow Fe Kq emission ( |Terashima et al.||2002[ ) and 
a handful of LINERs display broad dou ble-peaked Ha 
lines (e.g., Storchi-Bergmann et al.| [2003); these proper 



ties of the emission-line spectrum are consistent with the 



2 



Ncmmen et al. 



absence of a thin accretion disk, or a thin accreti on disk 
whose inner radius is truncated at > 100GM/c 2 ( Chen, 



Halpern fc Filippenko|1989||Chen fc Halpern|1989| . Last 
but not least, with the typical tuel supply of hot diffuse 
gas (via Bondi accretion) and cold dense gas (via stellar 
mass loss) available in nearby galaxies, LLAGNs would 
be expected to produce much higher luminosities than 
observed on the assumptio n of standa rd thin disks with 
a 10% radiative efficiency (Ho 2009). Taken together, 



ful outflows and jets, as suggeste d by theoretical stud- 



ies includin g analytical theory ( Narayan & Yi 1994 
M eier||20 0il |Nemmen et al.||2007[) and numerical sim 
ulations (McKimicy & Gammic 2004; Hawlcy & Kro- 
|hk| )2006[ fl'chekhovskoy, Narayan fc McKinney||2010[ ))' 
This is in line with the different observational studies 



which demonstrate t hat LLAGNs are generally radio- 



2008) 



loud |Ho 2002 200 8]) and accompanied by powerful jets 
'eg. 



Merloni & Heinz 



Hernz, Merloni fe Schwab| 2007| 
in fact, the so-called "radio mode" of AGN feed- 
back invoked in semi-analytic and hydrodynamic simula- 
tions of galaxy formation (e .g., Bowe r et al .|2006l|Croton 



et al. 



2006 Sijacki et al. 



this set of observational properties favors the scenario 
in which the accretion flow in LLAGNs is advection- 
dominatcd or radiatively inefficient. 

Advecti on-dominated accretion flows (ADAFfj for re- 
views see |Narayan Mahadevan fc Quataert 1998| |Yuan| would" presumably correspond to the AD Ah' accretion 

" state actively producing jets as explicitly incorporated 

in some wor ks ( Sijacki et al.|200"7 Okamoto, Nemmen & 
Bower|[2008 ). 



2007; Fanidakis et alT] 



)roton| 
|2011| ) 



2007[ | Narayan fcM cClintock 2008) are very hot, geo 
metrically thick, optically thin hows which are typified 
by low radiative efficiencies (L <C O.lAfc 2 ) and occur 
at low accretion rates (M < 0.01M E dd)- SMBHs are 
thought to spend most of their lives in the AD AF state 
( [Hopkins fc Hernquist|[2006l |Xu fc Cao|[2010|), the be st 



studied individual case being Sgr A* (e.g., Yuan||2007| ) 

In many LLAGNs, another component in the accretion 
flow besides the ADAF is required in order to account 



SEDs (lLasota et al. 1996 


Quataert et al. 1999 


Nem- 


men et aT. 2006 


Yu, Yuan & Ho|201l|) and as mentioned 



et al. 



Storch i-Bcrgm ann et al.|2003| |Eracleous, Lewis fc Flohic] 
2009[ ): the emission from a thin accretion disk whose in- 
ner radius is truncated at the outer radius of the ADAF. 
In many LLAGNs the accretion flow may begin as a stan- 
dard thin disk but somehow at a certain transition ra- 
dius the accretion flow abruptly switches from a cold to 
a hot ADAF mode. The details of how this transition 
might happen are still not well understood ( Manmoto 



20001 [Yuan fc Narayan|[2004{ |JNarayan fc McClln-, 
tock|2008*| ^ but it seems to be analogous to the transition 
between the different spectral states in black hole binary 



system s (Rcmi llard fc McClintock|2006 Done, Gierlihski| 
fc Kubota||2007[). 



Maoz ( 2007 ) challenged the scenario of the central en- 
gines of LLAGNs consisting of ADAFs and truncated 
thin disks. Maoz argues that LLAGNs in LINERs have 
UV/X-ray luminosity ratios similar on average to those 
of brighter Seyfert Is and based on that observation he 
posits that thin disks extending all the way down to the 
radius of the innermost stable circular orbits (ISCO) per- 
sist e ven for LLAGNs, des pite their smaller accretion 
rates. |Y u, Yuan feU o] (120111) showed that the SEDs com- 
piled by |Maoz|(|2007p are naturally fitted by ADAF mod- 
els with the addition of a truncated thin disk and also dis- 
cussed on theoretical grounds why the ADAF model has 
a superior explanation capability than the pure thin disk 
model. The ADAF model is the only model that can nat- 
urally account for the set of observational properties of 
LLAGNs within a self-contained theoretical framework. 
Hence, it is the physical scenario adopted in this work. 

ADAFs are relevant to the understanding of AGN feed- 
back since they are quite efficient at producing power- 

4 In this paper, we consider ADAFs and radiatively inefficient 
accretion flows (RIAFs) to be the same kind of accretio n flow so- 
lution. For a clarification regarding the terminology, see |Yu, Yuan| 
|fc Ho|2011| 



It is clear that an understanding of the physical na- 
ture of LLAGNs in LINERs will shed light on the nature 
of black hole accretion, outflows and consequently black 
hole growth and AGN feedback in present-day galaxies. 
One of the best ways of exploring the astrophysics of 
black hole accretion and outflows is by using multiwave- 
length observations of black hole systems and compar- 
ing the spectra predicted by specific models with the 
data. The goal of this work is therefore to probe the 
physics of accretion and ejection in the LLAGN popula- 
tion, by modeling their nuclear, broadband, radio to X- 
rays SEDs which provide constraints on physical models 
for the emission of the accretion flow and the jet. Fur- 
thermore, with a large enough sample of systems, we 
can derive from the fits to the data the parameters that 
characterize the central engines and build a census of the 
"astrophysical diet" of low-state AGNs. 

In this work, we will present how the currently fa- 
vored model for the central engines of weakly accret- 
ing black holes is able to account for the observed 
SEDs of LLAGNs. We will also present the typical ac- 
cretion/ejection - or feeding/feedback - properties of 
LLAGNs in LINERs as inferred from our detailed SED 
fits. With these parameters, we will then be able to 
draw a connection between the SMBHs in nearby active 
galaxies and the one in the Milky Way, Sgr A*. We will 
also present the inferred production site of the continuum 
emission and other properties inferred from the fits. 

The structure of this paper is as follows. In Section [2] 
we describe the sample of 21 LLAGNs in LINERs that 
we used, the data and classifications. Section [3] describes 
the physical model that we adopted in order to interpret 
and fit the SEDs and derive the central engine parame- 
ters. In Section[4]we describe the model fits to the SEDs, 
explaining in more detail the SED models for 5 objects 
which we consider the "reference" or fiducial LLAGN ex- 
amples. Section [5] presents the model fits to the other 16 
SEDs in our sample. In Section [6] we describe the distri- 
butions of central engine parameters resulting from our 
detailed model fits. We present the average SED result- 
ing from our fits - including the predicted emission in 
wavebands which can be observed with future facilities 
such as ALMA and the James Webb Space Telescope - 
and the inferred nature of the broadband emission in each 
wavelength in |7] We discuss in [8] which model is favored 
to explain our LLAGN sample, in particular which com- 
ponent in the flow dominates the X-ray emission. We 
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conclude by presenting a summary of our results in [9] 
and comparing the SMBH "diet" of massive black holes 
in LINER AGNs and Sgr A*. 

2. SAMPLE 



Heinz & Di Matteo||2003 ) 



Eracleous, Hwang fc Flohic] ( p010| (hereafter EHF10) 
compiled a sample of 35 SEDs of LLAGNs found in LIN- 
ERs which include high spatial resolution optical and UV 
observations with the Hubble Space Telescope (HST), as 
well as X-ray observations with Chandra. Most of the 
SEDs studied by EHF10 have also high-resolution ra- 
dio (observed with the Very Large Array - VLA - or 
VLBA/VLBI) and infrared observations available. In or- 
der to select the best SEDs to be modelled from the sam- 
ple of EHF10, we applied the following selection criteria: 
we selected only those objects for which their X-ray flux 
is not an upper limit. We also demand that there is a 
determination of the black hole mass. These selection cri- 
teria leave us with 24 LINERs. In three of these objects 
(NGC 404, NGC 5055 and NGC 6500) the emission is 
likely to be dominated by a stellar population (EHF10). 
Hence we discard these three LLAGNs and we are left 
with 21 objects which are listed in Table [T] This table 
also lists their corresponding Hubble and LINER Types, 
black hole masses, bolometric luminosities and Edding- 
ton ratios. 

According to the quality of the available data, we can 
further classify our sample of LINERs in two groups. 
Group A contains the objects that have the data with 
the best quality available. For instance, the radio, op- 
tical and UV bands have each at least one data point 
which is not an upper limit to the luminosity. Obviously, 
the objects in this group provide better constraints to 
the accretion-jet models. Group B comprises the sources 
that provide not as good contraints to the models, be- 
cause they lack observational constraints in the radio to 
UV region of the spectrum. Some objects in this group 
have only upper limits to the flux in the UV, others have 
no data points or only upper limits in the radio band. All 
the objects in both groups have observations in the near- 
IR band which are treated as upper limits. We take this 
approach because the observations in the near-IR were 
taken with lower spatial resolutions, i.e. larger apertures 
that include considerable contamination from the emis- 
sion of the host galaxy. Therefore, near-IR observations 
can only be considered as upper limits to the nuclear 
emission. 

The masses of the SMBHs were estimated from the 
stella r velocity dispersions using the Mbh — cr relation 



ship (Ferrares e fc MerrittpOOOj |Gebhardt et al.||2000 



Tremame~et al.||20U2[ ), with the exception of JNGC 266 



whose black hole mass was estimated from the measured 
5 GHz radio and t he 2-10 keV X-ray luminositie s using 
the correlation of Merloni, Heinz fc Di Matteo| (20031 
(the "fundamental plane of black hole activity ), and 
NGC 3031 (M81) and NGC 4486 (M87) whose black hole 
masses were estimated from the stellar and/or gas kine- 



matics ( Bower et al.|2000"j|Devereux et al.|2003 



et al.|2(J09| see Tableflj). As EHF10 note, for t 



Gebhardt 
re objects 

with multiple mass determinations using different meth- 
ods, the estimated masses are consistent with each other 
within a factor of 2. Due to the intrinsic scatter in the 
fundamental plane, the mass estimated for N GC 266 is 
subject to an uncertainty by a factor of ~ 10 (Merloni, 



Ihe optical- UV data of all objects were corrected for 
extinction as discussed in EHF10. In order to compute 
the bolometric luminosities from the SEDs, EHF10 used 
two methods. For the objects with the best sampled 
SEDs, they computed Lboi by integrating the SEDs di- 
rectly, ignoring upper limit data points. The objects 
with well-sampled SEDs are NGC 3031, NGC 3998, NGC 
4374, NGC 4486, NGC 4579 and NGC 4594. For these 
objects they assumed that pairs of neighboring points in 
the SEDs defined a power law, integrated each segment 
individually and summed the segments to obtain £t>oi- 
From this set of best sampled SEDs they estimated the 
average "bolometric correction" from the 2-10 keV lu- 
minosity to the bolometric luminosity (Lboi = 50Lx), 
which they used to obtain Lboi f° r the remaining ob- 
jects, for which the SEDs are not as well sampled. These 
bolometric luminosities are listed in Table [TJ 

3. MODELS FOR THE ACCRETION FLOW AND 

JET 

We fit the observed broadband SEDs of the LLAGNs in 
our sample using a model which consists of three compo- 
nents: an inner ADAF, an outer truncated thin accretion 
disk and a jet. The components of the model are illus- 
trated in Figure [T] We describe here the main features 
of this model. 



Jet 



Mass-loss via 
ADAF winds 



Thin disk 




Figure 1. Cartoon illustrating the model for the central engines 
of LLAGNs. It consists of three components: an inner ADAF, an 
outer truncated thin disk and a relativistic jet. 



3.1. ADAF component 

The inner part of the accretion flow is in the form of an 
ADAF which is a hot, geometrically thick, optically thin 
two-temperatur e accretion flow, which has low radiative 
efficiency (e.g., | Narayan, Maha devan fc Quataert 1998 



Kato, Fukue & Mincshige 1998). ADAFs are character 



ized by the presence of outflows or winds, which prevent 
a considerable fraction of the gas that is available at large 
radii from being accreted onto the black h ole. This has 
been suggested by numerical simulat i ons dStone et al 
19991 IHawley, Balbus fc Stone||2Tj0T| |igumenshchev et 
al.| |2003[ |De Villiers, Hawley fc Krolik||2003MMcKinney 
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Table 1 

Sample of galaxies and their basic properties 51 



Geilcixy 


Hubble 


Distance 13 




LINER 




Lu i 


-f'bol/iEdd 


SED 




Tvnp 
rype 


(Mpc) 


( A f~a tu /Mr^\ 
V ili BH/ J "Q ) 


Tvnp 


(gtct g t°\C 


( tzrtr G t\d 






type c 


NGC 266 


SB(rs)ab 


62.4 


W 


7.6 f 


LI 


7.4 X 10 40 


2.2 X 10 42 


4 


x 10- 4 


B 


NGC 1097 


SB(rl)b 


14.5 


( 1 \ 


8.1 


LI 


4.3 X 10 40 


8.5 x 10 41 


5 


x 10~ 5 


A 


NGC 1553 


SA(rl)0 


17.2 


J) 


7.9 


L2/T2 


1.3 x 10 40 


3.8 x 10 41 


4 


x 10~ 5 


B 


NGC 2681 


SBA(rs)0/a 


16.0 


\ 2 ) 


7.1 


LI 


6.1 x 10 38 


1.8 x 10 40 


1 


x 10~ 5 


B 


NGC 3031 (M81) 


SA(s)ab 


3.6 


(■V 


7.8 


S1.5/L1 


1.9 x 10 40 


2.1 x 10 41 


3 


x 10" 5 


A 


NGC 3169 


SA(s)a 


19.7 


(V 


7.8 


L2 


1.1 x 10 41 


3.3 x 10 42 


4 


x 10~ 4 


B 


NGC 3226 


E2 


21.9 


(2) 


8.1 


LI 


5.0 x 10 40 


1.5 x 10 42 


1 


x 10" 4 


B 


NGC 3379 (M105) 


El 


9.8 




8.2 


L2/T2 


1.7 x 10 37 


5.1 x 10 38 


3 


x 10 -8 


B 


NGC 3998 


SA(r)0 


13.1 


(2) 


8.9 


LI 


2.6 x 10 41 


1.4 x 10 43 


1 


x 10" 4 


A 


NGC 4143 


SAB(s)0 


14.8 


(2) 


8.3 


LI 


1.1 x 10 40 


3.2 x 10 41 


1 


x 10" 5 


A 


NGC 4261 


E2-3 


31.6 


(2) 


8.7 


L2 


1.0 x 10 41 


6.8 x 10 41 


1 


x 10" 5 


B 


NGC 4278 


El-2 


14.9 


(2) 


8.6 


LI 


9.1 x 10 39 


2.7 x 10 41 


5 


x 10~ 6 


A 


NGC 4374 (M84, 3C 272.1) 


El 


17.1 


(2) 


8.9 


L2 


3.5 x 10 39 


5.0 x 10 41 


5 


x 10~ 6 


A 


NGC 4457 


SAB(s)0/a 


10.7 


(1) 


6.9 


L2 


1.0 x 10 39 


3.0 x 10 40 


3 


x 10~ 5 


B 


NGC 4486 (M87, 3C 274) 


E0-1 


14.9 


(2) 


9.8 


L2 


1.6 x 10 40 


9.8 x 10 41 


7 


x 10~ 6 


A 


NGC 4494 


El-2 


15.8 


(2) 


7.6 


L2 


9.2 x 10 38 


2.8 x 10 40 


6 


x 10~ 6 


B 


NGC 4548 (M91) 


SBb(rs) 


15.0 


(3) 


7.6 


L2 


5.4 x 10 39 


1.6 x 10 41 


3 


x 10~ 5 


B 


NGC 4552 (M89) 


E 


14.3 


(2) 


8.2 


T2 


2.6 x 10 39 


7.8 x 10 40 


4 


x 10" 6 


A 


NGC 4579 (M58) 


SAB(rs)b 


21.0 


('1) 


7.8 


L2 


1.8 x 10 41 


1.0 x 10 42 


1 


x 10" 4 


A 


NGC 4594 (M104) 


SA(s)a 


9.1 


(2) 


8.5 


L2 


1.6 x 10 40 


4.8 x 10 41 


1 


x 10" 5 


A 


NGC 4736 (M94) 


(R)SA(r)ab 


4.8 


(2) 


7.1 


L2 


5.9 x 10 38 


1.8 x 10 40 


1 


x 10" 5 


A 



a The information in this table was obtained from EHF10, see text. 

k The number in parenthesis gives the source and method of the distance measurement, as follo ws: (1) From the ca talog of |Tully | l |1988| l, 
determined from a model for peculiar velocities and assuming Hq — 75 km s — 1 Mpc -1 ; (2) From|Tonry et al.|l|200 l[>, w ho used the surface 
brightness fluctation method. Followin g | Jensen et ah] J2Q03|, the distance modulus reported by^To n r^*etal.| l|200 f} wa s corrected by 
subtracting 0.16 mag; (3) From Frccdm anet**aT| " |T994j[20uT^ who used Cepheid variables; (4) From |(?ava^z7^taT*]^2009[ l who used the 
Tully-Fishcr method. 

c Lx is the X-ray luminosity in the 2 — 10 keV range. 

d The bolometric luminosities of NGC 1097, NGC 3031, NGC 3998, NGC 4374, NGC 4486, NGC 4579 and NGC 4594 were estimated by 
integrating the SEDs. For all other galaxies Lbol was determined by scaling Lx as described in [J2] 
e Classification of the SED depending on the quality of the data sampling (sec Section [2}. 

^ The mass estimate for NGC 266 is subject to a large uncertainty (~ 1 dcx) since it is based on using the fundamental plane of [Merloni, | 
|Heinz fc Pi Matteo| j2003( (see text). 



fc Gammie| |2004[) and analytical work ( Narayan fc Yi 
1QM| ^landlord fc Begelman|[l999l jNarayan et al 



Quataert & Gruzinov||2000ap . In addition, observational 



studies suggest that winds coming from the in ner parts 



Tombesi 



of accretion fl ows are ubiquitous in AGNs (e.g , 
et al.|2010a|b I and it is natural to take them into account 
when studying the LLAGN population. 

In order to take this mass-l oss into account, we follow 
Blandford & Begelman ( |1999[ ) and introduce the param- 
eter s by 



M = M„ 



R 

Row 



(1) 



to describe the radial variation of the accretion rate M out 
measured at the outer radius of the ADAF, i? ou t- Follow- 
ing the results of numerical simulations of the dynamics 
of ADAFs, in our models we allow s to vary over the 
range < s < 1, taking into account the current uncer- 
tainty in the determination of the amount of gas that is 
lost via winds in the accretion flow. 

The other parameters that describe the ADAF solu- 
tion are the black hole mass Al; the viscosity parameter 
a; the modified plasma f3 parameter, defined as the ratio 
between the gas and total pressures, (3 = P g /Ptou S, the 
fraction of energy dissipated via turbulence that directly 
heats electrons; and the adiabatic index 7. In our cal- 
culations, we adopt a = 0.3, /3 = 0.9 and 7 = 1.5. Tra- 



diti onal ADAF models a dopted 8 to be small (8 < 0.01; 
e.g., |Narayan fc Yi| l995). On the other hand, it has been 
argued that the value of 8 can be potentially increased 
due to different physical processes - such as magnetic re- 
connection - that affect the heating of protons and elec- 



trons in hot plasma s (e.g., Quataert & Gruzinov 1999 
Sharma et al.| |2007|) . Given the theoretical uncert amty 



related to the value of 8, we allow it to vary over the 
range 0.01 < 8 < 0.5. 

The case in which 8 = 0.5 and s = 0.3 deserves spe- 
cial attention, since these values were obtained from fit- 



ting current ADAF models to the SED of Sgr A* ( Yuan. 
Quataert fc Narayan| |2003| Yuan, Shen fc Zhi- Qiang 
2006[ ). We note that while [Yuan, Yu fc Ho| ( |2009[ ) fix 
the values of 8 and s in their hts, we leave these parame- 
ters free taking into account the theoretical uncertainties. 
There are degeneracies b etween s and 8 when fitting the 
data with ADAF models (Quataert et al. 1999), such that 
different combinations of the values of these parameters 
may result in similar SEDs. 

The cooling mechanisms incorporated in the calcula- 
tions are synchrotron emission, bremsstrahlung and in- 
verse Comptonization of the seed photons produced by 
the first two radiative processes. Given the values of 
the parameters of the ADAF, in order to compute its 
spectrum we first numerically solve for the global struc- 
ture and dynamics of the flow, as outlined in | Yuan et al.| 
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( |2000[ ); |Yuan Quataert fe Narayanj ( |2003[ ). Obtaining 
the global solution of the differential equations for the 
structure of the accretion flow is a two-point boundary 
value problem. This problem is solved numerically using 
the shooting method, by varying the eigenvalue j (the 
specific angular momentum of the flow at the horizon) 
until the sonic point condition at the sonic radius R s is 
satisfied, in addition to the outer boundary conditions. 

There are three outer boundary conditions that the 
ADAF solution must satisfy, specified in terms of the 
three variables of the problem: the ion temperature Tj, 
the electron temperature T e and the radial velocit y v (or 
equivalently the angu lar velocity f2). Following Yuan, 
Ma & Narayan (20081, when the outer boundary of the 
ADAF is at the radius -R ut = 10 4 i?s (where Rg is 
the Schwarzschild radius) we adopt the outer bound- 



We adopt a jet model based on the internal shock sce- 



ary conditions T Q 



0.2T v ; r , T 0] 



0.19T vir and 



A ut = 0.2, where the virial temperature is given by 
T vir = 3.6 x 10 12 (R s /R) K, A = v/c s is the Mach 
number and c s is the adiabatic sound speed. When 
the outer boundary is at R out ~ W 2 Rs we adopt the 
boundary conditions T out ,i = 0.6T V i r , T out! c = 0.08T V j r 
and Aout = 0.5. After the global solution is calculated, 
the spectrum of the accretion flow is obtained (see e.g., 
Yuan, Quataert & Narayan 2003 for more details). We 
verified that if these boundary conditions are varied by 
a factor of a few, the resulting spectrum does not change 
much. 

3.2. Thin disk component 

Outside the ADAF there is an outer thin accretion 
disk with an inner radius truncated at i? t r = -Rout and 
extending up to 10 5 i?s- Therefore, in practice the outer 
radius of the ADAF corresponds to the transition radius 
to the thin disk. The other parameters that describe 
the thin disk solution are the inclination angle i, the 
black hole mass and the accretion rate M out (the same as 
the accretion rate at the outer boundary of the ADAF). 
In the cases where i? ut ~ 10 4 i?s we simply ignore the 
contribution of the thin disk spectrum. 

The thin disk emits locally as a blackbody and we take 
into account the reprocessing of the X-ray radiation from 
the ADAF. This reprocessing effect has only a little im- 
pact on the spectrum of the thin disk though, with the 
resulting SED bein g almost identical to that o f a stan- 
dard thin disk (e.g., Frank, King & Rainc 2002 1. 



3.3. Jet component 



The SEDs of LLAGN s a re generally radio- loud ( Ho 
19991 |Ho fc Peng|[200T] |Ho||2002[ |Terashima fe Wilson 



20031 EHFlOTTTnit see |Maoz||2007[ |Sikora, Stawarz fe 



asota||2007[ ). The radio prominence of these SEDs is 
usually explained by invoking the synchrotron emission 
of relativistic jets, since the the accretion flow does not 



served radio luminosity (e.g., Quataert et al. 


1999 Ul- 


vestad & Ho 2001; Wu & Cao 2005; Nemmcn ct al. 2006 


Wu, Yuan & Cao|2007; Yuan, Yu & Ho 2009). 


Some au- 



explained by the jet component (e .g., Falcke, Kording & 



Markoff|2004l|Markoff et al.|2008| ). We therefore include 



in our modelling the contribution from the emission from 
a relativistic jet. 



glows (e.g., 


iPiran 1999 


Spada et al. 2001 


& Narayan 2005 for more details). This r 



sec 



Yuan, Cui 



'his model has been 
adopted in previous works to underst and the broadband 



SEDs of X-ray binaries and AGNs QSpada et al.|[200T| 
Yuan, Cui fc Narayan| [2005| |Nemmen et al.p006| \Wa] 



Yuan & Cao||2007p. Accord ng to this jet model, some 
fraction of the material in the innermost regions of the 
accretion flow is transferred to the jet producing an out- 
flow rate Mj et and a standing shock wave in the region 
of the jet closest to the black hole is formed. This shock 
wave is created from the bending of the supersonic ac- 
cretion flow near the black hole in the direction of the 
jet. We calculate the shock jump (Rankine-Hugoniot) 
conditions to find the electron and ion temperatures of 
the plasma (T e and Tj). We find that the jet spectrum 
is not very sensitive to changes in T e and Tj, since the 
emission is completely dominated by non-thermal elec- 
trons (see below). We therefore adopt T.- L = 6.3 x 10 11 K 
and T e = 10 9 K in our calculations of the jet emission. 

The jet is modelled as having a conical geometry with 
half-opening angle <f) and a bulk Lorentz factor Tj which 
are independent of the distance from the black hole. The 
jet is along the axis of the ADAF and makes an angle i 
with the line of sight. The internal shocks in the jet are 
presumably created by the collisions of shells of plasma 
with different velocities. These shocks accelerate a small 
fraction £ e of the electrons into a power-law energy dis- 
tribution with index p. The energy density of accelerated 
electrons and the amplified magnetic field are described 
by two free p arameters, e e and eg. F ollowing |Nemmen] 



et al 
calcu 



( |2006[ ); |Wu, Yuan fe Cab| ( |2007[ ), we adopt in our 
ations of the jet emission the values cfi = 0.1 ra- 



dians, £ e = 10% and Tj = 2.3, which corresponds to 
v/c sa 0.9 (except in the case of M87, which has an in- 
dependent estimate of Tj available). Unless otherwise 
noted, we adopt i = 30°. Therefore, there are four free 
parameters in the jet model: Mj 0t , p, e e and eb- In our 
calculations we consider the synchrotron emission of the 
jet, with the optically thick part of the synchrotron spec- 
trum contributing mainly in the radio and the optically 
thin part contributing mainly in the X-rays. 

3.4. Free parameters and fit procedure 

When fitting the observed SEDs with our coupled 
accretion-jet model, we have eight free parameters. Four 
of these parameters describe the emission of the accretion 
flow: the accretion rate M out , the transition radius be- 
tween the inner ADAF and the outer truncated thin disk 
Rtr, the fraction of viscously dissipated energy that di- 
rectly heats the electrons S and the strength of the wind 
from the ADAF s. The other four parameters character- 
ize the jet emission: the mass- loss rate in the jet Mj e t, 
the electron energy spectral index p, and the electron and 
magnetic energy parameters e e and €b- M out and M- ]ct 
are obviously expected to be correlated, but given our 
ignorance about the mechanism of jet formation we vary 
these parameters independently in our fits. We require 
for consistency that Mj et /M out < 1. Since the ADAF 
loses mass via winds (equation [T]), it is also of interest to 

check the value of Mj ct /M(3Rs), in which the accretion 



G 



Nemmen et al. 



rate is calculated at the radius of the innermost stable 
circular orbit for a Schwarzschild black hole. 

Throughout this paper we will use the dimension- 
less mass accretion rates m = M/M^dd, noting that 
the Eddington accretion rate is defined as A/Edd = 
22M/(1O 9 M ) M yr _1 assuming a 10% radiative effi- 
ciency. We will also express the black hole mass in terms 
of the mass of the sun, m = M/M@, and the radius in 
terms of the Schwarzschild radius, r = R/Rg. 

Our fitting procedure can be summarized as follows. 
We vary the free jet parameters in order to fit the radio 
observations, since the radio is emitted predominantly by 
the jet. Once we fit the radio band with the jet model, 
we try to fit the truncated thin disk model to the IR- 
optical data, in order to constrain the transition radius 
and Af out (keeping in mind that the theoretical spectrum 
should not exceed the observed luminosities). We then 
use the X-ray data to constrain the ADAF component 
in order to refine the estimates of i?tr and Af out , and 
to estimate s and 8. Depending on the combination of 
parameters that we use, the jet can have an important 
contribution to the emission not only in the radio, but 
also in the X-ray band. We sum the emission of the 
inner ADAF, outer thin disk and jet, and compare this 
sum with the observed SEDs. 

We explore the parameter space of the models consid- 
ering the many possible fits to the data. We search the 
literature for independent constraints from other works 
on the values of the model parameters, such as the mass 
accretion rate, inclination angle and the transition ra- 
dius, as well as the jet power which is obtained from the 
jet model. 

4. MODEL FITS TO WELL SAMPLED SEDS 

We describe in this section the results obtained from 
fitting our coupled accretion-jet model to the SEDs that 
were selected using our selection criteria in Sect ion [2] W e 
describe in more detail the resul ts fo r NGC 4594 H j4.ll) , 
NGC 4374 {UM, NGC 4486 ({jO), NGC 3031 (p~4| 
and NGC 399S(f45j). These five objects will serve as 
the "reference" examples that will illustrate the general 
results for the whole sample. These objects were cho- 
sen for a more detailed discussion because of their good 
multiwavclcngth sampling. We show the model fits to 
the other 16 SEDs and the resulting model parameters 
in Section [3 

In the SED plots that follow below, the error bars rep- 
resent the uncertainty in the extinction corrections ap- 
plied to the the optical-UV data by EHF10. The lower 
bars correspond to the measurements without any ex- 
tinction correction, while the upper bars represent the 
same observations after a maximal extinction correction 
(see EHF10). The arrows represent upper limits to the 
nuclear luminosity. These upper limits are caused by 
either non-detections, or because the corresponding ob- 
servations were taken at lower spatial resolutions. In 
the latter case, the upper limits include a potentially 
significant contamination by the emission from the host 
galaxy. This happens particularly for the observations in 
the near-IR band. 

Table [2] lists the model parameters that result from 
the SED fits for the AD and JD type of models. We 
list whenever available the Bondi accretion rate and the 



corresponding reference from which meondi was taken. 
Jj°j s corresponds to the jet power estimated either from 
observations or using t he correlation between je t power 

pmod 
-Tjet 



Merloni & Heinz 



(2007) 



and radio luminosity of 

represents the jet power resulting from our jet model 



calculated as P 



jet 



4.1. The case of NGC 4594 



Pellegrini (12005 1 studied Chandra X-ray observations 
oi JNGC 4594 and derived the Bondi accretion rate 
msondi = 2±?-| x 10- 3 . 

The SED of NGC 4594 is plotted in Figure M together 
with two different spectral fits. The left panelof Figure 
[2] shows an ADAF-jet model in which the ADAF com- 
ponent (dashed line) completely dominates the emission 
from the IR to the X-ray bands and is consistent with the 
available optical-UV and X-ray data. In this model, the 
jet (dot-dashed line) dominates the radio emission and 
contributes only weakly to the X-ray flux. The solid line 
corresponds to the sum of the ADAF and jet emission. 

The parameters of the ADAF model are rh out = 
6.3 x 10 -3 (consistent with the Bondi accretion rate es- 
timated by |Pellegrin7||2005| ), r out = 10 4 , 8 = 0.01 and 
s = 0.3. The jet parameters are rhj e t — 9 x 10~ 7 , p = 2.3, 
e e = 5 x 10~ 3 and eg = 0.03. The jet kinetic power is 
Pjet/Lboi ~ 4. In other words, the kinetic power carried 
by particles in the jet exceeds the overall radiative power. 
The ratio of the mass-loss rate in the jet to the accretion 
rate estimated at 3R$ is mj et /m(3i?s) rj 1.6 x 10~ 3 . In 
other words, only ss 0.2% of the mass that is ultimately 
accreted by the black hole is channeled into the jet. 

The dotted line in the left panel of Figure 12] shows an 
ADAF + truncated thin disk model for whicn r tr = 300, 
m out = 9.1 x 10~ 4 , 8 = 0.01 and s = 0.1. The thin 
disk emission corresponds to the "bump" at 10/im. This 
model demonstrates that accretion flows with a much 
smaller transition radius are also consistent with the 
data. Since we have only upper limit points in the mid 
to near-IR for NGC 4594, the available data are not suf- 
ficient to constrain the presence of a truncated thin disk 
in this source. 

The right panel of Fig. [2] shows a model for the SED 
of NGC 4594 in which the jet dominates the radio and 
X-ray emission. In fact, the right panel shows that the 
relativistic jet can account well for the entire observed 
SED. The jet parameters are m- }e t — 4.5 x 10 -7 , p — 2.01, 
e e = 0.8 and e B = 3x 10~ 2 , with P iet = 8.4xl0 42 <**" 



erg s 



Pjet/iboi « 4 and m- iCt /m{m s ) « 0.009. 

We display in the right panel of Figure [2] one illustra- 
tive ADAF model which is consistent with the possibility 
that the jet emission dominates in this source, i.e. the 
ADAF contribution in all bands (except at 1 mm) is very 
small compared to the jet. The ADAF parameters are 
?Tiout = rhsondi) r out = 10 4 , 8 = 0.01 and s = 0.3. We can 
see that the ADAF contribution in this case is only sig- 
nificant in the wavelength ranges ~ 1 mm — 100 /im. The 
importance of the ADAF emission in the other bands is 
small compared to the jet. 

By exploring the parameter space of the accretion/jet 
models within the range of plausible values allowed by 
theory, we demonstrate that there are basically two pos- 
sible types of models which can accommodate the ob- 
served SED of NGC 4594 and the other LINERs in 
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Figure 2. Models for the SED of NGC 4594 (Sombrero). The dashed and dot-dashed lines show respectively the emission from the 
ADAF and jet, while the solid line shows the sum of the radiation from these components. Left: model in which the ADAF dominates 
the observed X-ray emission. The dotted line shows for comparison an ADAF + thin disk model in which the transition radius is 300-Rs. 
Right: model in which the jet dominates the X-ray output. See the text for the parameters. 



this section. In the first type, the emission from the 
ADAF dominates the observed X-rays; in the second type 
of model, the jet emission dominates the X-rays. We 
will hereafter use the abbreviations AD (as in ADAF- 
dominated) and JD (as in jet- dominated) when referring 
to the former and latter types of models, respectively. A 
third type of model is also possible in which the jet and 
the ADAF contribute with similar intensities to the high 
energy emission. These results apply not only to the five 
LINERs discussed in more detail in this section but also 
to the other sources in our sample with a few exceptions. 

4.2. NGC 4374 

A prominent jet resolved with the VLA is observed 
to create cavities in the X-ray emitting gas (Allen et al. 
2006 Finoguenov et al.||2008|). The Bondi accretion rat e 



is m Bo ndi = 4 x 10~ 4 "( |Pellegrini|2005| | Allen et al.||2006| ). 

The left panel of Figure [5] shows an ADAF-jet model 
(dashed line) in which the ADAF completely dominates 
the emission from the near-IR to the X-ray bands and 
explains well the available optical-UV and X-ray data, 
as was the case for NGC 4594. The synchrotron peak 
of the ADAF also dominates the emission in the range 
~ 1 mm — 100 /im. The jet (dot-dashed line) dominates 
the radio emission, but contributes weakly to the X-ray 
flux. The thin disk dominates the IR emission and its 
contribution peaks at A ~ 10 /im, creating a "red bump" 
in the SED as op posed to the "big blue bum p" observed 



the jet power estimated by |Merloni k, Heinz (2007) based 
on the calorimetry of the X-ray cavities observed with 
Chandra (3.9 x 10 42 erg s _1 ). The ratio rhj e t/m(3Rs) is 
1.5 x 10~ 3 . 

The right panel of Fig [3] shows a model in which the 
jet dominates the radio and the 1-100 keV X-ray band. 
The emission in the region between 1000 A and 1 keV is 
dominated by the Compton peak of the ADAF spectrum. 
As was the case for the previous model, the IR emission 
is dominated by a truncated thin disk with a small tran- 
sition radius (rt r = 30) which is even smaller than the 
AD fit. At v > 100 keV the bremmstrahlung emission of 
the ADAF is the dominant radiation process. The pa- 
rameters of ADAF are rh ou t = 1.5 x 10~ 4 (in agreement 
with the Bondi rate) and S = s = 0.3. The jet parame- 
ters are Th\ e t — 1.6 X 10 -6 , p — 2.2, e e = 9 x 10~ 3 and 

3.4 x 10 42 erg s" 1 « 17L bo i- 



cb 



'jet - 

8 x to" 3 , 



with P\ 



The ratio rhj e t/rh(3Rs) is rh- ]e t / 'rh(3Rg) ~ 0.02. There 
are no data available to constrain the predicted ADAF 
bumps at soft and hard X-rays in this model. Notice 
that, similarly to the JD fit to the SED of NGC 4594, the 
shape of X-ray spectrum from the ADAF with S = 0.3 
(dashed line) does not agree with the data. 

4.3. NGC 4486 



in quasars (e.g., Koratkar & Blaes 1999 Nemmen & 



Brotherton||2010 ) . The parameters of the accretion fio 



are m ou t = 3.95 x 10" 4 , r tr = 150, S = 0.01 and s = 0.1. 
As was the case for NGC 4594, the mid to near-IR data 
are insufficient to make the case of a truncated thin disk 
compelling and similar models with much larger values 
of rt r are not ruled out by the available IR data. 



The jet parameters are m 



jet 



io- 



P 



2.4, 



Di Matteo et al. ( 2003 1 previously calculated the Bondi 
accretion rate using Chandra X-ray data and found it 
to be M B o„di ~ 0.1 M G yr" 1 (m Bc ,„di « 7 x lO" 4 ). 
They also fitted the mult i wavelength SED of M87 with 
an ADAF model for different values of S (but not in- 
cluding mass-loss in the accretion flow, i.e. s = 0) 
and f ound that the ADAF emission rough ly explains the 
SED. |Biretta, Sparks fc Macchetto| (p99) analysed HST 
observations of the jet in 1V187 and estimated tha t T,- > 6 
and 10° < i < 19°. Base d on the results of |Biretta^ 



0.01 and €b = 0.1. The resulting jet power is Pj C t = 2.1 x 



Sparks & Macchetto 
elling the parameters 



1999), we adopt in our jet mod- 



6 and i = 10°. The kinetic 



10 erg s ~ 4Lboi which is in excellent agreement with power carried by the jet is estimated to be in the range 
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Figure 3. Models for the SED of NGC 4374/M84 showing the emission of the ADAF (dashed), jet (dot-dashed), truncated thin disk 
(dotted) and the total emission (solid). Left: model in which the ADAF dominates the observed X-ray emission ("AD model"). Right: 
model in which the jet dominates the X-ray output ("JD model"). 



1Q43 _ 1Q 44 g -i ( Bicknell fc Bcgclman 1999 
Owen, Eilek fc Kassim||200(j| |Allen et al .||lW0(j| |Merloni 



39; 



& Heinz 2007), therefore the jet power resulting from 
fitting our jet model must be within this range. 

The left panel of Figure [4] shows an ADAF-jet model 
in which the ADAF dominates the optical-UV (OUV) 
and X-ray emission (i.e. AD-type model), while the jet 
dominates the radio band. The parameters of the ac- 
cretion flow are m out = 5.5 x 10~ 4 (in agreement with 
the Bondi rate), r out = 10 4 , 5 = 0.01 and s = 0.1. 
The jet parameters are rhj e t = 5 x 10~ 8 , p = 2.6, 
10~ 3 and €b — 8 x 10 -3 . The jet power is given 



Yu&Ho 
JD mode 



( 2009 ) instead successfully fit the data with a 
using p = 0.5. 

4.4. NGC 3031 

This source is among the brightest LLAGNs known 
since it is the nearest AGN besides Centaurus A and 
has been the subject of a broadband mu l tiwavelength 



moni toring campaign (Markoff et al. 
20101 



Due to the 



2008 Miller et al. 



by 



jet 



= 6 X 



with Pjct/iboi ~ 6 and 



e B = 

10 42 erg s" 

m jet /m(3i? s ) ~ 9 x 10~ 5 . 

The right panel in Fig [4] shows a jet model which ex- 
plains quite well the radio and X-ray observations (JD- 
type model), although it underpredicts the OUV data 
by a factor of a few. The jet parameters are mj c t = 
3 x KT 8 , v = 2.3, e P = 10~ 3 and e R = 8 x 10" 3 , with 



many available observational con- 
straints, accretion-jet models can be well constrained for 
this LLAGN. For illustrative purposes, here we will re- 
port our modelling results using the data compiled by 
EHF10. 



P jet = 9 x 10 42 erg s" 1 and P ]ot /L 
The SED of M8 7 was previously fitted using ADAF 



bol " 



and /or jet models 


Di Matteo et al. 


2003 Yuan, Yu&Ho 


2009 

TTTTTV 


Li et al. 2009 


. Di Matteo et a 

' — A 1 "i A 1 \ II — 


. (|2003) modelled the 



Devereux fc Shearer] ( 2007 ) modeled the profile of the 
broad double-peaked Ha line with a relativistic thin disk 
model. They were able to explain the profile with an 
inclination angle of 50° . Such a high inclin ation angle is 
supported by radio obs ervations of the jet ( |Bieten holtz, 
Bartel fc Rupen||2000[ ). The Ha line profile of M81 is 
consistent with an inner radius of ~ 280 — 360i?,c; fo r the 
line-emitting thin disk (Devereux & Shearer 2007). In 



of 5 but not including mass- loss (i.e. 
w,r IT ; ni l donnnh 



, s = 0). The model 
adopted by |Li et al. l ( I 2009I ) is quite similar to that of 
Di Matteo et al.] ( |2003p althoug h the former inco r porat e 
general relativ istic corrections. Di Matteo et al. (2003); 
Li et al. ( 2009[ ) obtained that the ADAF emission with no 
mass-loss approximately reproduces the SED and results 
in accretion rates consistent wi th the Bondi rate. Our 
AD model is similar to that of Di Matteo et al. (20031 
but it also incorporates the effect of winds, as suggested 
b y numerical simu l ations 



Yuan Yu fc Ho| ( |2009[ ) tried to model the SED using 
an ADAF model with d = 0.5 but failed to fit to data 
with an AD model. The reason is that for high values of 
S (S > 0.1) and small accretion rates (r7i out <§C 0. 01) the 
ADAF X-ray spectrum is harder than the data. |Yuan,| 



our models for the SED, we therefore adopt an inclination 
angle i = 50° and we impose the constraint on the models 
that r tr < 360. 

The left panel of Figure [5] shows the AD fit model for 
the SED of M81. In this model, the synchrotron radi- 
ation from the jet is the dominant mechanism of radio 
emission. In the wavelength range between « 10 pan and 

1 fim a "red bump" can be seen which corresponds to 
the emission of the truncated thin disk. The thin disk 
dominates the emission in this band. In the interval be- 
tween 1 mm and « 100 /urn the ADAF is the dominant 
component. The ADAF is also the flow component that 
dominates the emission for wavelengths less than 1 /i un- 
til X-rays. The sum of the radiation of all components is 
plotted as the solid line. The accretion flow parameters 
are m out = 3 x 10 -3 , r tr ~ 360 (consistent with the Ha 
modeling results), S = 0.01 and s = 0.16. The jet param- 
eters are mj ct = 2x 10~ 6 , p — 2.2, e e — 0.1 and es — 0.01. 
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Figure 4. Same as Figure[3]for M87/NGC 4486. 
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Figure 5. Same as Figure |3]for M81. 



The jet power is given by Pj 



jet 



x 10 41 erg s 1 with 



1(T 



-Pjet/£boi ~ 4 and m jet /m(3i?s) 

In this AD model, a low value of S = 0.01 is required in 
order to reproduce the X-ray spectral shape character- 
ized by r = 1.88. As was the case for M87, if we increase 
the value to 8 = 0.3 keeping the other parameters fixed 
except s — 0.4 the predicted value of the photon index is 
reduced and the ADAF spectrum becomes harder than 
the data. 

The right panel of Fig. [5] shows a JD model in 
which the jet component dominates the radio, optical- 
UV and X-ray bands of the SED. In the wavelength 
range between 1 mm and ss 100 /im the ADAF is the 
dominant component, although there are no observa- 
tional constraints to the model in this region. Between 
« 10 fim and 6000 A the truncated thin disk dominates 
the emission as is the case for the AD model. The pa- 
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rameters of the accretion flow are m out = 8 x 10 , 
rtr = 50, 5 = 0.3 and s = 0.6. The parameters of 
the jet are m- ]ct = 1.2 x 10~ 5 , p — 2.05, e e = 0.6 and 
e B = 10- 4 (P je t = 4.8 x 10 42 erg s" 1 , P jct /L hol a 23 and 
m ict /m{3Rs) « 3 x 10" 5 ). 

We adopted two different transition radii in the fits 
shown in figure [5] The accretion flow model with r tr = 
360 overpredicts the UV data (left panel of fig. [5 1. A 
smaller transition radius of r tr = 50 on the othernand 
does a better overall job of describing the optical-UV 
data. Such small transition radius seems also to be fa- 
vored by the modelling of the Fe Ka emission line ( |Young| 
|2007[ ). We note however that with the available 
R and 



et al. 



near- 



optical data it is not possible to set a strong 
constraint on the properties of the truncated thin disk. 

M81 has been the subject of several modeling efforts 
dQuataert et al.||1999l iMarkoff et al. 1120081 [Yuan, Yu & 
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Ho|2009 ) . Our goal here is not to discuss extensiv ely this 
particular LINER since this was done elsewhere (Miller 
et al. 2010), but rather put the general features of our 
AD and JD models for M81 in context with the results 
f ound by other authors 



UV data. The parameters of the jet are m 



jet 



3.5 x 



Quataert et al. ( 1999) did not consider the jet contribu- 
tion and favored an AD scenario for the origin of the SED 
(they used an "old" A DAF model assuming S = 0.01 and 



0, i.e. no winds). |Markoff et al.| ( |2008| ) ; |Yuan, Yu fe 
Ho| ( 2009 ) on the other hand were also able to expl ain the 
observations, bu t in the context of a JD model. |Yuan, 
Yu & Ho ( 2009 ) also considered the simultaneous contri- 
bution of the ADAF and jet but restricted t hemselves to 
5 = 0.5 in their fits. Mainly for this reason, Yuan, Yu & 
Ho| ( |2009[ ) were unable to find an AD model for the SED, 
since as we mentioned before, ADAF X-ray spectra with 
such values of 6 are quite h ard. 

It should be noted that |Quataert et al. ( 1999| ) esti- 
mated a transition radius which is twice (rt r = 100) the 
value favored in the present work and a much higher 
accretion rate. The reason for this is that these authors 
previously fitted the SED of M81 using an outdated value 
of the black hole mass, which was fifteen times smaller 
and more uncertain than the value we use (estimated via 
spatially resolved gas and stellar kinematics, see EHF10). 
As a consequence, the multicolor blackbody radiated by 
their thin disk for the same accretion rate and transition 
radius that we use is considerably hotter and fainter. As 
a consequence, they needed to increase r tr and m. 



4.5. NGC 

The shortest wavelength (1750 A) UV data point in 
the SED of NGC 3998 presented in EHF10 is ar toma- 
lously high because of variability (Devereux 2011). As 
discussed by Devereux this data point was obtained many 
years before all the other observations, when the source 
was much brighter. For this reason, we exclude this UV 
measurement from our analysis and plots. 

The left panel of Figure [6] shows an AD model in which 
the ADAF dominates the emission from the IR to the 
X-ray bands and explains well the available optical-UV 
and X-ray data. As usual, the jet dominates the radio 
emission. The ADAF parameters are m ou t = 7.2 x 10~ 3 , 
r out = 10 4 , S = 0.1 and s = 0.4. The jet parameters are 
m jet = 1.7 x 10" 6 , p = 2.2, e e = 0.01 and e B = 10~ 3 . 
For the resulting jet power from this model (Pj e t = 9 x 
10 42 erg s" 1 , m ]et /rh(3Rs) « 6 x 10~ 3 ), the jet kinetic 
energy is quite modest compared to the radiative output 

(Pjet/ibol « 0.6). 

In the AD model above we adopted a large outer radius 
for the ADAF, but smaller radii are not ruled out by the 
data. For instance, we are able to obtain a reasonable 
fit to the SED with r tI = 500, m out = 10~ 3 and s = 
S = 0.01. This model is consistent with the IR upper 
limits and accounts for the X-ray data. The transition 
radius cannot be much smaller than this value, otherwise 
the emission of the truncated thin disk woul d exceed the 
IR up per limits. Furthermore, as noted by |Ptak et al 



(2004), the lack of Fe Ka line emission also suggests that 
the value of r tr is not so small. 

The right panel of Fig. J6] shows a JD model in which 
the jet dominates the radio, optical-UV and the X-ray 
emission. Note that although the jet accounts quite well 
for the X-rays, it somewhat underestimates the optical- 



10 



P 



2.01, 



1.8 x 10 43 



e e = 0.75 and e B = 3 x 10~ 5 (P jot = 
erg s -1 , Pjot/iboi ~ 1-3 and m jQt /m(3Rs) « 
0.025). For illustration, we also show in the right panel 
an ADAF model computed with parameters such that 
its contribution in X-rays is weak (m ou t = 4 x 10~ 4 , 
r tr = 100, S = 0.3 and s = 0.3). Note that, as for 
the cases discussed in the previous section, for values of 
S > 0.3 the ADAF X-ray spectrum is harder than the 
data. 

The SEP of th is NGC 3998 was modelled before by 
Ptak et al. ( 2004 1 . Here we confirm the main results of 
their work, namely that both AD and JD models are able 
to account well for the broadband SED. 

5. FITS TO SPARSELY SAMPLED SEDS 

We present in this section the model fits to the 16 SEDs 
that were not discussed in the previous sect ion. 

In the sub- 



The SED fits are shown in Figures 7]|T2 



sections below, we briefly describe the aetails of the fits 
to the SED of each object. The model parameters are 
summarized in Table [2] 

5.1. NGC 0266 

The AD model accounts slightly better for the best-fit 
X-ray slope, but note the pronnounced uncertainty on 
the value of observed photon index. The JD model re- 
quires p < 2 which is below the usual r ange 2 < p < 3 
suggested by relativistic shock theory ( Bednarz & Os- 
trowskifl998| |Kirk et aLpOOO] ) . 



5.2. NGC 1097 

Note that even though the SED of NGC 1097 is as 
well sampled as many of the LLAGNs in Section |4j we 
chose to include it in this section rather than discussing 
it in more detail in j|4] because its SED was already well- 
studied by |Nemmen~et aL] ( |2006[ ). 

We take into account here the HST UV spectrum which 
was not included in EHF10 because it contains contri- 
butions from an obscured nuclear starburst in addition 
to the AGN. Thi s spect rum was obtained from |Storchi-| 
Bergmann et al. ([2005) who studied and modeled it in 
detail. 

Our models do not take into account the contribution 
of the nuclear starburst. Therefore, the AD and JD mod- 
els do not reproduce the "UV bump" between s» 1000 A 



and ~ 7000 A (see Fig. 3 in Storchi-Bergmann et al 



2005) 



The AD model displ ayed in Figure [7 corre sponds to 
the model obtained by Nemmen et al.j l2006). The JD 
model is also able to account for the whole SED, though 
it does not reproduce the slope of the X-ray spectrum 
as well as the AD model. Regarding the JD fit, given 
the available data there is no need to incorporate the 
contribution of the ADAF. Hence we do not include any 
ADAF model in the right panel of Fig. [7j 

5.3. NGC 1553 

The Bondi accretion rate was estimated by | Pellegrini] 
( [2005] , 

Even though the JD model roughly accounts for the es- 
timated Lx , it fails to fit the slope of the X-ray spectrum 
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Figure 6. Same as Figure|3]for NGC 3998. 

even with a small value of p. We take this as evidence 
that this source is unlikely to be JD. The estimated m out 
is consistent wit h the lower limit on WBondi obtained by 
Pellegrinil (12005 1. 



5.4. NGC 2681 

For this LLAGN, there are not enough optical obser- 
vations to fit the truncated thin disk model and estimate 
the transition radius. 

5.5. NGC 3169 

As was the case of NGC 2681, for this LLAGN there are 
not enough optical observations to fit the truncated thin 
disk model and estimate the transition radius. The JD 
model accounts well for the available data. The inferred 
extreme values of e e and e b imply that essentially all the 
energy in the post-shock region of the jet is carried by 
the particles. 

5.6. NGC 3226 

The SED of this LLAGN is well fitted by both an AD 
and JD type of models. As is the case of NGC 2681, there 
are no good optical band constraints on the emission of 
the truncated thin disk. 

5.7. NGC 3319 

This source has only two data points outside the X-ray 
band, one in the radio and the other in the optical, both 
of which correspond to upper limits. Therefore, there are 
few constraints for th e accretion-jet mod el. The Bondi 
rate was estimated by David et al. (2005). 

Figure [9] shows that the AD model is able to reproduce 
the observed SED given the few constraints available, 
but there are no radio data to constrain the jet model 
in this case so we don't include it. The available data 
are not enough to constrain well the transition radius. 
The accretion rate required by the AD model is more 
than an order of magnitude higher than meondi- One 
possible explanation for this result is that the accretion 
rate is enhanced by gas released by stars. This is in line 
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with the findings of Soria et al. ( 2006a|b I for a sample of 
quiescent early- type galaxies. 

5.8. NGC 4143 

The contribution of the truncated thin disk is required 
in order to account for the optical emission. The required 
transition radius in particular is r tr = 70. 

The JD model requires a value of p smaller than the 
usual values suggested by relativistic shock theory, as was 
the case of NGC 266. 

5.9. NGC 4261 

For this object, estimates of the Bondi rate, inclination 
angle and jet power are available dGliozzi, Sambruna fc 
Brandt|[2003) |Merloni fe Heinz||2007p . There are two in- 
dependent estimates ot the jet power. The first is based 
on energetic considerations and the VLBA observations 
of the pc-scale jet, which yield a lower lim it to the jet 
power (Gliozzi, Sambruna & Brandt 2003). In the sec- 
ond i jne_tEoT]wi|lis^ and 
the Merloni & Heinz ( 2007 ) correlation obtaining a much 



higher jet power estimate. 

Note that the accretion rate that we obtain from the 
AD fit is somewhat smaller than msondi- This AGN is 
likely to be strongly affected by extinction in the OUV 
band (EHF10). Hence, the OUV measurements probably 
do not capture the emission of the central engine. For 
this reason, it is not surprising that the models overpre- 
dict the OUV emission by almost one order of magnitude 
in Fig. [9| 

5.10. NGC 4278 

The radio emission of NGC 4278 w as previously stud- 
ied in the c ontex t of ADAF models by Di Mat teo, Carilli 



& Fabian ( |2001[) , who found that s > is required in 



order tor the radio emission predicte d by ADAF models 
not to overestim ate the radio flux. Di Mat teo, Carilli| 
& Fabian (20011 also es timated that the Bon di rate is 
CTUOr- 0.01. iGiroletti et all (|2005l found a 



m B 



two-sided radio structure for the jet using VLA data 
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Figure 7. SEDs and coupled accretion-jet models for NGC 0266, NGC 1097 and NGC 1553. The left panel shows the AD models for 
each object while the right panel displays the JD models. The dashed, dotted and dot-dashed lines correspond to the emission from the 
ADAF, truncated thin disk and jet, respectively. The solid line when present represents the sum of the emission from all components. 
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Figure 8. Same as Figure[7]for NGC 2681, NGC 3169 and NGC 3226. 
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Figure 9. Same as Figure[7]for NGC 3379, NGC 4143 and NGC 4261. 
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and estimated that the jet is oriented close to the line of 
sight (2° < i < 4°) and mildly relativistic (T 3 ~ 1.5). In 
our jet models we therefore adopt i — 3° and Tj — 1.5. 

Both the AD and JD models require a truncated thin 
disk with a small transition radius (r tr ~ 30 — 40) in 
order to account for the near-IR data. The values of the 
transition radius that we used in the models presented 
above (r tr = 30 — 40) are not unique and models with 
larger radii can also reproduce the SED. For instance, 
an accretion model with r tr = 100, m out = 4 x 10~ 3 , 
S = 0.1 and s — 0.77 is able to explain the optical and 
X-ray observations. 

The jet emission by itself underpredicts the 1 mm radio 
observation, therefore the ADAF contribution is required 
even in the JD case. In the JD fit, a small value of s is 
favored by the observations as is the case of the JD fit 
for NGC 4143. 

5.11. NGC 4457 

Since there are only upper limits in the radio band, 
the jet power for this object was estimate d using the 
observation at v = 1.5 x 10 10 Hz and the Merloni & 



emission. The truncated thin disk in our fits is some- 
wh at hott er than the thin disk modeled by |Quataert et 



Heinz (2007) correlation. 



5.12. NGC 4494 
Same as NGC 4457. 

5.13. NGC 4548 

Same as NGC 4457. Given the large uncertainty in the 
photon index, the X-ray spectrum does not provide good 
constraints for the models. 

5.14. NGC 4552 



The Bondi accretion rate was estimat ed by |Merloni fc 



Heinz| ( |2007) (see also |Allen et al.||2006| from the X-ray 
profiles ot density and temperature! The kinetic power 
carried by the j et was estimated by |Merloni fc Heinz] 
(2007) (see also Allen et al. 20061 from the energy de- 
posited in the X-ray cavities. 

Since the X-ray spectrum is quite soft, AD models are 
unable to account for the X-ray emission. The JD model 
in Fig. ITT] is roughly consistent with the radio observa- 
tions and explains quite well the X-ray data, but overpre- 
dicts the optical data. The resulting j et power is roughl y 
consistent with the valu e estimated by Allen et al. (2006) ; 
Merloni & Heinzl 120071). 



5.15. NGC 4519 

NGC 4579 shares many characteristics with NGC 3031 
(M81) and NGC 1097: its nucl eus features broad d ouble- 
peaked Balmer emission lines ( |Barth eta l. 2001) and a 
lack of the iron Ka line emission ( |Eracleous et al.| [2002 ) . 
Quataert et al. (19991) previousl y modeled the SEP of 
NGC 4 579 cc^iIeTT^ [Ho| ([ T999l) using the "old" ADAF 



model ( jNarayan fc Yi|r9"9"S b. 
From the width of the broad Ha line, IBarth et al. 



(2001) obtained a rough estimate of the inner radius of 



the line- emit ting portion of the accretion disk r tr ~ 160. 
Hence, we adopt in our SED models r tr < 160 and 
i = 45°. 

As can be seen in Fig. |11| the optical continuum data 
do not require a small transition radius since it can be 
accounted by either the jet (JD fit) or ADAF (AD fit) 



al. (1999), since their adopted black hole mass is out- 
dated and is ~ 10 smaller than the value we use in this 
work. The transition radius can be chosen to be some- 
what smaller if the accretion rate is lower and still be 
consistent with the data. 

We show together with the JD model an ADAF model 
which is the same as the one calculated for the AD fit 
except for the higher value of s. 

5.16. NGC 4736 

As was the case of NGC 1097, NGC 4143 and NGC 
4278, this LINER requires a small transition radius in 
order to explain the OUV data. 

6. THE DISTRIBUTION OF THE ACCRETION- JET 
PARAMETERS FOR LLAGNS IN LINERS 

Our detailed SED modeling allows us to place indepen- 
dent constraints on the mass accretion rate onto the black 
hole, the jet mass-loss rate and jet kinetic power for the 
LLAGNs in our sample. These quantities are of direct 
relevance to studies of the feeding and feedback of the 
underfed supermassive black holes in nearby galaxies. It 
is therefore prudent to study the resulting distributions 
of the central engine parameters from our models. 

6.1. Accretion rates 

In the previous section we showed that essentially two 
types of models (AD and JD) are able to account for the 
SEDs, depending on the combination of parameters that 
we choose. We found that in most of the JD fits, the con- 
tribution of an underlying ADAF is not required at all 
(cf. the JD models for NGC 4594, M87 and NGC 4579). 
Since the estimate of the accretion rate relies on fitting 
the ADAF-thin disk model to the data, in those cases no 
robust to constraint is available. We therefore show in 
the left and right panels of Figure |T3| the histograms dis- 
playing the distribution of values oTm ou t (dimensionless) 
and Mout (in units of M Q year -1 ) respectively, consider- 
ing only the AD models. For one object (NGC 4552) we 
are unable to find a suitable AD model, therefore we do 
not include this object in the histogram. 

Most of the objects have accretion rates in the range 
m « 0.003 — 0.01. The mean accretion rate for the sam- 
ple, in the case of the AD fits, is (log to) = —2.5 or 

log ^M/M Q year -1 ^ = —2 with standard deviations 

of 0.5 and 0.6, respectively. 

Not all gas supplied at the outer radius of the accre- 
tion flow ends up being accreted onto the black hole due 
to mass-loss via winds produced in the ADAF. From 
our derived parameters we can estimate the fraction of 
the gas that actually falls into the event horizon. Fig- 



14 



mg 



shows the histogram of values of M(r; sco ) tak- 
mto account the radial dependence of the density 
profile in the ADAF, where r isco = 3. We obtain 

(log (A/ (r isco )/M year 1 ) ) = -3.3 with a standard 

deviation of 0.5. In other words, typically only a frac- 
tion of 14% of the available fuel supply reaches near the 
ISCO and presumably a similar fraction is accreted onto 
the black hole. 
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Figure 10. Same as Figure [7] for NGC 4278, NGC 4457 and NGC 4494. 
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Figure 11. Same as Figure [7] for NGC 4548, NGC 4552 and NGC 4579. 
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Figure 12. Same as Figure [7] for NGC 4736. 
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Figure 13. Distribution of accretion rates at the outer radius of the ADAF obtained from our AD SED fits, a, left: Accretion rate in 
Eddington units, b, right: Accretion rate in physical units. 



6.2. Jet powers 

Especially relevant for studies of the feedback of mas- 
sive black holes in the local universe is the amount of 
energy that the jets deposit in their environment. Fig- 
ure [15] shows the distribution of jet kinetic powers ob- 
tainedvia both the AD and JD type of SED fits. As we 
noted above, we were not able to estimate a jet power for 
the AD model fit to NGC 4552. The average jet power 
log (Pjet/erg s -1 )) is 41.8 and 42.2 respectively for the 
D and JD scenarios, each with a standard deviation of 
ss 1 dex. Therefore, we can say that the typical jet power 
is ~ 10 42±1 erg s _1 . The distributions in Fig. [l5]are in 
qualitative ag reement with the distribu t ion of jet powers 
estimated by Nagar, Falcke & Wilson ( 2005[ ) (see their 
Fig. 6). 

Also relevant in studies of the kinetic AGN feedback 
is what fraction of the rest mass energy associated with 
the mass accretion rate is tapped by the central engine 
and converted into the jet kinetic power, which can be 



defined as 77j Gt = P- iC t/{Mc 2 ). The choice of accretion 
rate in the definition of the "jet kinetic efficiency" is ar- 
bitrary and we could choose to use either the accretion 
rate M ut at the outer radius of the ADAF or the one at 
the innermost stable circular orbit (ISCO), Mi 8CO (note 
that Mi sco < M ou t given the mass-loss through winds in 
the ADAF). The ISCO radius in this case corresponds to 
3i?5, appropriate for a Schwarzschild black hole. Figure 
[T6| shows the distribution of the jet kinetic efficiency de- 
fined using -Mout (top panel) and Mj sco (bottom panel) 
for the case of AD SED fits, since in the JD fits there 
are no reliable estimates of the accretion rates. The av- 
erage values of log J7j c t are -2.9 (M = M ont ) and -2.1 
(M = M isco ) (both these averages have an uncertainty 
of Ri 0.7 dex). Clearly, the latter efficiency is higher in 
order to be consistent with the smaller ISCO accretion 
rates. 

We are also able to quantify the fraction of the mass 
accretion rate that escapes the attraction of the black 
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Figure 14. Distribution of accretion rates at the radius corre- 
sponding to the ISCO of a Schwarzschild black hole, estimated 
from the AD fits. 
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Figure 16. Distribution of jet kinetic efficiencies Pj ot /(Mc 2 ) con- 
sidering M = Mj sco and M = M ou t for the AD models. 




Figure 15. Distribution of jet powers. 

hole and is channeled into the jet. Figure [17] shows 
the distribution of values of M- iCt / Mi BCO for the AD 
fits. As it turns out, the average value of the ratio is 

log ^Mjct/ 'M\scoj^ = —2.9, i.e. on average 0.4% of the 

material that reaches the ISCO is channeled in the jet. 
The corresponding uncertainty in this estimate is 0.6 dex. 

6.3. Bolometric luminosities and radiative efficiencies 

One possible way of estimating the bolometric lumi- 
nosities for the LLAGNs using the observed SEDs is to 
integrate the SEDs directly, not taking into account up- 
per limits (especially in the IR) and assuming a suitable 
interpolation. For instance, this procedure was carried 
out by EHF10 who also computed a bolometric correc- 
tion for the X-ray luminosity. 

In our case, we have physically motivated fits for 
the 21 SEDs. Therefore, for each LINER in our sam- 
ple, we calculated Lboi by integrating the total SED. 
Carrying this out using the SEDs obtained according 
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Figure 17. Distribution of Mj ct /Mi B 



to the AD scenario we obtain the average value of 
(log (^boi/ erg s -1 )) = 41.5 with a standard deviation 
of 1 dex. For the SED fits computed according to 
the JD scenario we find a corresponding average of 41.3 
with a similar standard deviation, which gives the same 
average as the one computed using the bo lometric lu- 
minosities estimated by EHF10. Figure 18 displays the 
distribution of values of Lboi resulting from the models 
compared with the one estimated by EHF10. 

Defining the X-ray bolometric correction as kx = 
Lboi/Lx, we calculated the average bolometric correc- 
tion (logKx) resulting from our models as 1.45 and 1.35 
according to the AD and JD scenarios respectively, with 
standard deviations of 0.36 dex (AD) and 0.46 dex (JD). 
The corresponding geometric means of Kx are 28 (AD) 
and 22 (JD), the mean values are 43 (AD) and 47 (JD) 
and the medians are 22 (AD) and 14 (JD). Our models 
result in values of kj which are somewhat lower than 
those estimated by EHF10. 
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Figure 18. Distribution of bolomctric luminosities inferred from 
the models compared with the ones estimated by EHF10. 




Figure 19. Distribution of the ratio of jet power to bolomctric 
luminosity. 



Using our SED models we can calculate what is the 
typical fraction of the accreted matter that is converted 
into the combined radiation of the accretion flow and 
the jet, i.e. the radiative efficiency r) Ta< i = L\ m \/ (M out c 2 ) 
(note that the accretion rate in the definition of this ef- 
ficiency is the one at the outer radius of the ADAF). 
We obtained the average value (log^ ra d) = —3.25 with a 
standard deviation of 0.75 dex. Hence, the typical radia- 
tive efficiency is 5 x 10~ 4 . 

We estimated how the jet kinetic power compares to 
the total radiative output of the central engine, i.e. we 
calculate the ratio Pj e t/-^boi where Pj et is obtained from 
the jet model and Lboi from integrating the modeled 
SEDs. We obtain (log (P jet /L bo i)) ps 0.2 for the AD-type 
of fits and 0.9 for the JD models (there is a ps 0.5 dex 
uncertainty in both values). The reason for the smaller 
jet powers in the AD fits compared to the JD ones is 
that the mass-loss rates required in the former fits are 
systematically smaller than in the ones in which the jet 
has a higher degree of contribution in the X-rays. Based 
on our results and given the degeneracy in the fits, we 
can say that the typical kinetic to radiative power ratio 
of LLAGNs is ps 5 (i.e. the average of the value from 
the A D and JD fits). This is in agreem ent with the find- 
ings of Nagar, Falcke & Wilson (2005) that the accretion 
power outp ut in LLAGJNs is dominated by the jet power. 
Figure 19 shows the distribution of values of Pj e t/Pboi 
resulting from the models. 

7. THE AVERAGE SED 

It is of interest to compute the average SEDs resulting 
separately from the AD and JD scenarios. These SEDs 
are useful for different purposes. First and foremost, 
since the observed SEDs are the main observables that 
we use to derive the central engine parameters via our 
fitting, as a consistency check the average model SEDs 
should reproduce the average SED obtained directly from 
the observed ones. 

Secondly, SEDs obtained from averaging observed data 
points are obviously limited by the observed bands. Out- 
side the observed bands different authors adopt ad-hoc 



interpolations between t he data, u sually a linear interpo- 
lation in log-log space |Hoj|1999] EHF10), which might 
not reflect the actual physical processes involved in the 
emission. The average from the AD and JD scenarios 
that we obtained can hence serve as guidelines - with a 
physical justification - to the typical shape of the SEDs 
in the bands which have not been constrained yet. 

Motivated by the reasons above, we show in Fig. [20^ , 
the average SEDs computed separately for the AD and 
JD scenarios. We first normalized the individual SEDs 
to the same X-ray luminosity of 10 40 erg s _1 in the band 
2-10 keV. This value is approximately the average X-ray 
luminosity for the LINERs in our sample. After normal- 
izing the SEDs, for each frequency bin we computed the 
geometric mean of the model SEDs corresponding to the 
particular model class as (log vL v ). Following EHF10 we 
choose to compute the average of the logarithm of the lu- 
minosities instead of using the values of vh v themselves 
in order to reduce the effect of outliers in the resulting 
average. 

The different bumps in the average model SEDs reflect 
the different physical processes that operate in the flows. 
In the average AD model SED, the different bands are 
dominated by different radiative processes in the flow as 
described below: 

1 m 1 mm: jet optically thick synchrotron emission. 

1 mm — ps lOO^im: the first bump and usually the 
strongest one in the average AD SED is due to the 
ADAF synchrotron emission. 

ps 10/im — ps 1/im: the second bump is due to the ther- 
mal emission from the truncated thin disk. 

ps 1/Ltm — 0.1 keV: there is a weak UV bump due to 
the first inverse Compton scattering of synchrotron 
photons in the ADAF. 

ps 0.1 keV — 100 keV: the last bump in X-rays is dom- 
inated by the second inverse Compton scattering of 
synchrotron photons in the ADAF. 
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The average JD SED is more simple and has only one 
bump in the IR between a few xlO/zm - a few xlfim 
which is due to the truncated thin disk. The rest of 
the emission is due to synchrotron photons from the jet. 
Overall, the spectral shapes of the average JD and AD 
SEDs arc quite similar in the radio and X-rays above 
1 keV. In between these bands the shapes are slightly 
dissimilar. The bumps in the average SEDs are much 
less pronounced than the bumps in the individual SEDs. 
This occurs because the bumps in the individual SEDs 
do not peak at the same place and when the SEDs are 
averaged these bumps are smoothed out. For this reason, 
the average model SEDs will not resemble any one of the 
individual SEDs. 



Figure 20 1 also displays the average data points com- 
puted from the observed SEDs in a similar way by EHF10 
where the error bars represent the uncertainty in the 
emission due to the uncertainty affecting the amount 
of extinction correction involved. The shaded region 
around the average best-fit X-ray power-law represents 
the standard deviation in the value of the LINER pho- 
ton index. In order to illustrate the wide diversity of the 
individual SEDs, we show in Figure [20fc the la scatter 
affecting the model SEDs (where we show only the AD 
average SED for simplicity) and the observed ones. 

As expected, the average model SEDs agree well with 
the observed constraints. There are some details that 
are worth mentioning. The shape of the X-ray spectrum 
of the average JD SED is slightly softer than the corre- 
sponding shape of the average AD SED. Both are within 
the la uncertainty in the photon index of the average ob- 
served SED. In the OUV, even though the model SEDs 
agree with the observed constraints they are quite dif- 
ferent from each other. For instance, the red bump is 
stronger in the average JD SED. The average JD SED 
predicts a lower level of UV flux. In the radio band, the 
model S EDs are quite similar to each other. 

Figure 21 shows the average AD and JD SEDs com- 
pared to the average ones of radio-loud a nd radio-quiet 
quasars computed by iShang et al.| (|20ilj). The average 



quasar SEDs computed by Sha ng et al.| (2011 ) are very 
similar to the ones of |Elvis| ( |l994| but the former in- 
clude more detailed features, and are based on more re- 
cent da ta o btained with improved instrumentation. As 
in Fig. 20 1, the SEDs were normalized such that they 
all have the same X-ray luminosity in the 2-10 keV band 
of 10 40 erg s" 1 . The UV excess in the quasar SEDs (the 
big blue bump) is clearly apparent in comparison to the 
LLAGN ones. It is also interesting that for v > 10 17 Hz 
the average AD, average JD and the radio- loud quasar 
SEDs are quite similar. 

8. DISCUSSION: THE NATURE OF X-RAY 
EMISSION 

In our AD models, the X-rays are produced predom- 
inantly by inverse Compton scattering by the ADAF of 
seed synchrotron photons produced in the accretion flow 
itself; in contrast, in our JD models, the X-rays are dom- 
inated by the optically thin synchrotron emission in the 
base of the jet. In this regard, the Chandra X-ray spec- 
trum is one of the most important constraints for the 
origin of high-energy photons in LLAGNs. 

There are only two objects for which the AD or the 
JD can be ruled out, within the observational uncertain- 



ti es a nd the explored mo del pa rameter space: NGC 1553 
( p!3| and NGC 4552 (j pTl4| ). For NGC 1553, the JD 
model does not reproduce the hardness of the X-ray spec- 
trum whereas for NGC 4552 the AD model is not able 
to account for the corresponding X-ray softness. For the 
remaining objects, by inspecting the SED fits we can see 
that both the AD and the JD models are able to repro- 
duce with varying degrees of success the whole SED from 
radio to X-rays, in particular the X-ray portion. 

In order to quantify which model systematically fits 
better the broadband SEDs we use as a goodness-of-fit 
parameter a modified \ 2 . We define this modified y 2 
simply as 



X 



(yob 



2/model J 



(2) 



where y = log vL v and the sum is carried over the ob- 
served data points. We sample the X-ray continuum us- 
ing 3 data points from 0.5 keV to 10 keV and we do not 
take into account upper limit data points. The values of 
X 2 for each model are displayed in Table [2j 

Considering the 20 objects for which we were able to 
obtain both JD and AD fits, we find that \ 2 {JD) < 
X 2 {AD) for 9 objects and x 2 (JX>) > X 2 {AD) for the 
remaining 8. In other words, for 9 objects (10 including 
NGC 4552 for which we did not find an appropriate AD 
fit) the JD results in better SED fits as opposed to 8 for 
which the AD model is a better fit. The results from 
this goodness-of-fit test show that 52% of the SEDs in 
our sample are better fitted by the JD model whereas 
38% are better reproduced by the AD scenario. For the 
3 objects with X 2 (JD)/ X 2 (AD) « 1 (NGC 1097, NGC 
3998 and NGC 4374) it is hard to tell - based on the 
broadband SED alone - which model provides a better 
fit. 

We should point out a number of caveats involved when 
considering which scenario better accounts for the nature 
of the X-ray emission. From an observational point of 
view, even though we are carrying the largest systematic 
analysis of the SEDs of LLAGNs in LINERs done up to 
date, the size of our sample is limited. We need to extend 
our modeling to a larger number of LLAGNs in order to 
increase the statistics and draw further conclusions. Sec- 
ondly, for many sources we only have a few data points 
available to fit, so clearly for these objects we need to 
obtain more measurements to better constrain the mod- 
els. In addition, in our modeling we find that one of 
the main observables that help to pinpoint the dominant 
source of X-rays is the hardness of the X-ray spectrum. 
ADAF models tend to predict harder X-ray spectra while 
jet models predict softer ones. In most of the LLAGNs 
in our sample there is a significant observational uncer- 
tainty in the X-ray hardness derived from the Chandra 
observations which make it difficult to determine which 
whether the AD or JD is favored based on the X-ray 
spectrum. 

From a theoretical perspective, there are considerable 
theoretical uncertainties involved in the the ADAF-jet 
models. For instance, the values of the model param- 
eters that affect the microphysics and dynamics of the 
ADAF and jet are not sufficiently constrained. In par- 
ticular, 6 and p have a large impact on the hardness of 
the X-ray spectrum predicted by ADAF/jet models but 
there is a considerable uncertainty in the range of pos- 



22 



Nemmen et al. 



44 
43 
42 
41 
40 
39 
38 
37 
36 
35, 



lm 10cm 1cm 1mm 100)im 10cm 1/ttn 1000.4 .IkeV IkeV lOkeVlOOkeV 



AD mean 
JD mean 



— ■ 



10 



12 



14 

log(i//Hz) 



16 



18 



20 





44r 




43- 




42- 


7 


41; 






cn 


40- 


a) 




""- 




jjf 


39 - 


Dl 




O 


38; 




37- 




36- 











lm 10cm 1cm 1mm lOO^m 10/m 1/jm lOOOA.lkeV IkeV lOkeVlOOkeV 




10 



12 



14 

logW Hz) 



16 



18 



20 



Figure 20. a, left: The average SEDs (geometric mean) computed separately for the AD and JD models (dashed and dotted lines 
respectively). The data points correspond to the geometric mean computed by EHF10. b, right: 1<t scatter around the average (model 
and observed) SEDs illustrating the diversity of individual SEDs. The solid line shows the average AD SED and the shaded region 
corresponds to the standard deviation from the AD models. The points correspond to the mean computed by EHF10 and the error bars 
show the scatter in the measurements. In the OUV, the filled circles correspond to measurements without any reddening correction whereas 
the open circles correspond to the maximal extinction correction. 
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The nature of the X-ray emission in LLAGNs has been 
debated in the last few years by several authors favoring 
in some cases the AD or J D models based on the analysis 



of individual sources (e.g., Falcke fc Markoff 2000 Yuan 



Markoff & Falcke 2002 Yuan et al. 2002; Yuan. Quataert 



fc Narayan||2003| |Wu, Yuan fc Cao||2007t |Markorf et al 



2008| |Miller et aL} 2010) or in a statistical sense based 
on the fundamental piai 



plane of black hole activity flMer 
loni, Heinz fe Pi Matteo|2003| |Falcke K5rding fc Markoff 
20041 [Yuan fc Cui||2005| [Yuan, Yu fc Ho||2009[ |Plotkin 
et al.||201ip. We suggest different routes to clarify this 



issue in the future. 

The AD and JD models should predict different c harac- 
teristic rad io and X-ray variability timescales (e.g.JPtak 
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Figure 21. The average JD and AD model SEDs compared to 
the average radio-l oud and radio-quiet quasar SEDs computed by 
|Shang et af] | |2011[ l. 

sible values of these parameters despite r ecent progress 



et al.fl998 ) , hence we should be able to test which com 
ponent is dominant in X-rays by carrying out a simul- 
taneous monitoring of the variability of radio and X- 
rays. By comparing the variability pattern predicted by 
jet/ADAF models we could pinpoint the X-ray dominant 
component. This promising strategy is quite similar to 
what the observational campaigns carri ed out for M81 
( |Markoff et~aT1|2008[ |Miller et aLl|2010D a nd the radio- 
quiet Seyfert 1 NGC 4051 ( [King et al.|20ll[ ) and could be 
applied to many othe r LLA gJNs. We note that in the case 



(e.g., Sharma et aL]|2007 Medvedev 



20061). Considering 
the overall goodness-ot-tit of the AD and JD models, it 
is important to point out that our jet model is basically 
a phenomenological one. Uncertainties affecting for ex- 
ample the jet model could translate in larger values of 
X 2 - For instance, when we find \ 2 (JD) < x 2 (AD), that 
could be possibly because the JD model is less well con- 
strained than the AD one, i.e, we have more freedom in 
the jet fitting compared to the ADAF model. 

Taking into account these caveats, it is still premature 
to favor one model over the other based only on the SED. 



of M81, |Miller et al 
X-rays originate in t 



(2010) favor a scenario in which the 
le transition region between the thin 
NGC 4051 is a different case since 

0.05L Edd 



disk and the ADAF; 

it is a relatively bright AGN accreting at the ■ 
level (Peterson et al. 2004) and its jet-disk coupling is 
quite distinct compared to LLAGNs. 

The Fermi Gamma-ray Space Telescope can be quite 
helpful in this regard . LLAGNs are potential sources of 
7-rays (Takami 20111 and in fact one of the source in 
this sample (NG C 4486) has alrea dy been detected by 
the Fermi LAT ( |Abdo et aL]|2009[ ). By using the model 
parameters that we derived in our radio-to-X-rays fits we 
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should be able to predict the 7-ray spectr um (Mahade- 
van, Narayan fe Krolik|1997 Takami|2011 ) and compare 
with tuture Fermi detections. In this way we should be 
able to compare AD/JD predictions and through 7-ray 
observations place further constraints on the production 
site responsible for the high-energy emission and the jet- 
disk connection in LLAGNs. 

Observations in the mm and sub-mm are also very 
helpful because they constrain the ADAF synchrotron 
emission and hence also the X-ray emission. This fol- 
lows because the synchrotron photons in the ADAF are 
inverse-Compton-scattered to X-rays. 

Finally, further refinements and a better unders tand- 
ing of the fundamental plane of black hole activity (Mer 



and JD scenarios is ~ 10 42 erg s 1 ranging from 



10 40 to 



10 44 erg s" 1 . 



loni, Heinz fe Pi Matteo|2003||Falcke, K ording & Markoft 
2004|]Gultekin et al.|2009| [Yuan, Yu fc Ho|2009| |Plotkin 
et al.||201ip should make the fundamental plane a better 
tool for constraining the radiative processes shaping the 
radio and X-ray emission in sub-Eddington black hole 
source and by consequence LLAGNs. 

9. SUMMARY 

We performed detailed modeling of the broadband - 
radio-to-X-rays - spectral energy distributions of a sam- 
ple of 21 low-luminosity AGNs in LINERs selected from 
EHF10. With this exploratory modeling, our goal is to 
constrain the general properties of the central engines of 
the dominant population of AGNs at z « 0. Our coupled 
accretion-jet model consists of an accretion flow which is 
radiatively inefficient in the inner parts and becomes a 
thin disk outside a certain transition radius. The rela- 
tivistic jet is modeled in the framework of the internal 
shock scenario. 

We demonstrated that there are two classes of mod- 
els that can explain the majority of the observed SEDs. 
We call the first one AD which stands for "ADAF- 
dominated" since the ADAF dominates most of the 
broadband emission, particularly the X-rays. In the sec- 
ond class of models, the jet component dominates the 
majority of the continuum emission and for this reason 
we call this scenario JD as in "jet-dominated" . The SEDs 
predicted by the AD models have a more complex shape 
than the JD ones, given the richer variety of radiative 
processes involved in ADAFs. For both scenarios though, 
the radio band is almost always dominated by the syn- 
chrotron emission from the jet and the near IR to optical 
band is dominated by the truncated thin disk thermal 
emission. 

From our exploratory modeling of the SEDs we are able 
to constrain important parameters that characterize the 
central engines of LLAGNs, as summarized below. 

Mass accretion rates: based on the AD models, the 
values of the accretion rates supplied to the accre- 
tion flow lie in the range 4 x 10~ 4 — 0.02 in Ed- 
dington units with a mean (logm) = —2.5 or alter- 
natively (log (m/Mq year" 1 ^ = -2 (standard 

deviations of « 0.5 dex respectively). Of that gas 
supply, typically only ~ 10% is accreted onto the 
black hole with the remaining gas being lost due to 
outflows. 

Jet powers: the typical jet power resulting from the AD 



Jet production efficiencies: The typical efficiency 
with which the rest mass energy associated with 
gas supplied to the accretion flow is converted into 
jet kinetic power is estimated to be ~ 0.001 (range 
~ 10~ 4 — 0.01). The efficiency is ten times larger if 
we consider the gas that reaches a few gravitational 
radii only. 

Fraction of material ejected in the jet: On aver- 
age, a fraction of 0.4% of the gas that reaches a few 
gravitational radii is chanelled in the relativistic 
jet (range 0.01 - 1%). 

Radiative efficiencies: The typical efficiency of con- 
version of rest mass energy associated with gas sup- 
plied to the accretion flow into disk+jet radiation is 
«5x 10~ 4 (dispersion of 0.8 dex), producing an av- 
erage bolometric luminosity of a few xlO 41 erg s _1 . 

Kinetic to radiative output ratio: LLAGNs typi- 
cally release ~ 5 times more kinetic power in the 
jets than the total energy radiated away. 

Transition radii between thin disk and ADAF: A 

truncated thin disk is required to fit the available 
mid and near- IR data for only 4 so urces - NGC 
1097 (see also |Nemmen et~at||2006[ ), NGC 4143, 
NGC 4278 and NGC 4736 - with the resulting 
transition radii in the range 30 — 225Rs- 

The values that we derived from our SED models above 
provide useful indicators of the feeding and feedback 
properties of LLAGNs in LINERs, and by extension of 
the whole LLAGN population. 

Even though we used the physical scenario that is fa- 
vored to account for the combined set of observational 



prope rties of LLAGNs (|Ho T2008| |Narayan fe McClintock 
2008 ) , the inferred values of the parameters are model de- 
pendent. There are theoretical and observational sources 
of uncertainty in the values of the derived parameters. 
For instance, there are degeneracies between the effect of 
the model parameters on the resultin g SEDs ( Quataert & 
Narayan 1999b Nemmen et al. 2006 1. The uncertainty m 



the OU V - due to uncertain degree of extinction correc- 
tion - and the X-ray data also contribute to uncertainty 
in the estimated parameters. Furthermore, our sample is 
currently limited to 21 objects. In spite of these caveats, 
the estimates we obtained for the properties of central 
engines of LINER AGNs should be illustrative of the gen- 
eral behavior of the class. 

We can draw a link between the supermassive black 
holes in LINERs and the quiescent black hole in our 
galaxy, Sgr A*. By modelling the SED of Sgr A* with 
ADAF models, an accretion rate of m out ~ 10~ 5 was es 



timated (PYuan, Quataert fc Narayan 2003 Yuan, Shen & 
Zhi-Qiang 2006 ) . This value is two orders of magnitude 
below the typical accretion rate of LINERs that we esti- 
mated. Therefore, we could say that in order to Sgr A* 
"light up" and become a LINER-like LLAGN, it needs 
to accrete at a rate ~ 100 times higher than the present 
one. 
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The individual and average SED models can be useful 
for a different number of applications. They provide a li- 
brary of templates for the compact emission of LLAGNs, 
which can be useful for studies of the emission of other 
astrophysical systems in the nucle us, i.e., stellar popula- 
tio ns, dust and PAH feature s (e.g., Storchi-Bergmann et 



al. 2005 Mason et al. 2007). Furthermore, the models 



provide the expected emission on bands that were not 
yet observed in many LLAGNs (e.g., mid and far IR, far 
UV). 
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